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Kinematic outliers in the LMC: constraints on star-star 
microlensing 
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ABSTRACT 

Although a decade of microlensing searches towards the Large Magellanic Cloud 
(LMC) has detected 13 — 25 possible microlensing events, the nature and the location 
of the lenses, being either halo machos or LMC stars, remains a subject of debate. 
The star-star lensing models generically predict the existence of a small population 
(more than ~ 5%) of stars with a spatial and kinematic distribution different from 
the thin, young disc of the LMC. Here we present the results of a large spectroscopic 
survey of the LMC, consisting of more than 1300 radial velocities measured accurately 
with the 2dF instrument. In this large sample, no evidence is found for any extraneous 
population over the expected LMC and Galactic components. Any additional, kine- 
matically distint, population can only be present at less than the 1% level. We discuss 
the significance of this finding for the LMC self-lensing models. 

Key words: gravitational lensing - dark matter - galaxies: halos - galaxies: individual 
(LMC) - galaxies: kinematics and dynamics 



1 INTRODUCTION 

About a decade of searches for microlenses along the 
line of sight to the Magellanic Clouds (MCs) by MA- 
CHO/OGLE/EROS and a number of follow-up surveys have 
found more than 13 candidate microlensing events towards 
the LMC and two t owards the SMC. The MACHO collab- 
oration (Alcock et 



al 



200C) conclude that the lenses are 
mostly machos residing in the halo of the Galaxy, perhaps in 
the form of ancient white dwarfs (WDs). These WDs would 
make up 10-20% of the halo dynamical mass, with the re- 
maining 80-90% in some other comp onent. The detection 
of high proper motion white dwarfs (Hamblv^t^L 19971 



Star-star lensing models typically invoke an unvirialised 
component of the LMC and the SMC (e.g., a puffed up disc 
or a hot stellar halo with a surface brightness only a few 
percent of the LMC disc, a wrapped-around tidal ring, an 
offset bar etc, all due to tidal shocking among the LMC, 
SMC and the Galaxy) because of the low self-lensing optical 
depth of a thin virialised LMC disk (Gould 1995). Any un- 
virialised component of the LMC can then lense with stars 
i n the thin disc of the LMC. Some recent theoretical models 
(£hao & Evans- pOOlj |jetzer et al] |2002j, e.g.,) argue that 



Hodgkin et al. 



2000 



[bata et al 



200C; Oppenheimer et al 



2001) supports this interpretation. 

The assertion that WDs ma ke up at least 3% of t he ex - 



pected local dark matter density Oppenheimer et al 



Richei 



2001) 



( |2001| ) 

Some 



has been hotly contested (see review by 
argue very strongl y that these WD s trace the thick disk 
kinematics instead ( fr^eid et"aL 2001). If the thick disk inter- 
pretation is true, then the halo would have very few WDs, 
but the big question remains: what objects are responsible 
for microlensing the LMC stars? 

The answer is perhaps the LMC stars themselves. From 
the very beginning there have been many varia tions of star - 
star self-lensing proposals as first advocated by |3ahu| ( [l994| ). 



there could be enough stellar lenses in the LMC bar and 
disk to account for from half to all of the observed events. 
These models generally predict that the events should have 
peculiarities in photometry, kinematics and spatial distribu- 
tion because the lensing optical depth is higher in some re- 
gions than the others. These are highlighted in at least five 
of the dozen observed events: the LMC near-clump event 
MACHO-LMC-1 and LMC binary events MACHO-LMC- 
9 and MACHO-LMC-14, the SMC caustic event MACHO- 
98-SMC-l and long duration event MACHO-97-SMC-1. Al- 
though there are no strong direct evidences, it is reasonable 
to expect that all of the dozen observed lenses belong to the 
same LMC or SMC population as for these exotic events 
because ordinary single stars should be at least as common 
as binary stars in a given LMC or SMC stellar population, 
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and they all have similar lensing cross-sections and detection 
efficiencie s. 

Zhao flL999a| Jb|) discussed a number of observational 
tests to differentiate the two competing classes of lensing 
model. In particular, stars in any proposed unvirialised com- 
ponent should be identifiable as kinematic outliers, meaning 
that they deviate from the rotation curve of the LMC, as 
traced by the majority of LMC disc stars and by HI gas. 
The objective of the present contribution is to analyse the 
kinematics of large sample of stars towards the LMC to in- 
vestigate whether outliers can be singled out. These kine- 
matic outliers may reveal the presence of any hot compo- 
nent around the LMC, or any cold halo stream either in the 
Galactic halo or the LMC halo (Graff et al. 2000). We aim 
to be able to detect a kinematically distinct population (po- 
lar ring or thickened disk) in the LMC present at the ~ 5% 
level. If such a population is identified their relation to the 
microlensed sources can be studied. 

We foresaw two possible outcomes of this experiment: 

(i) The LMC microlensed sources and the randomly se- 
lected LMC field stars have disk kinematics: This would 
imply that the lenses are in the halo, supporting the MA- 
CHO interpretation of a dark baryonic component in the 
halo in the form of white dwarfs and consistent with the 
high- velocity, blue, old white dwarfs in the solar neighbour- 
hood detected by Ibata et al. (2000) and Oppenheimer et 
al. (2000). 

(ii) The sample of microlenses or the sample of field stars 
contain many outliers from the LMC rotation curve, imply- 
ing that the Galactic dark halo is almost entirely made of 
a non-compact component. The solar neighbourhood white 
dwarfs are merely a local peculiarity. 



2 OBSERVATIONS 

A sample of ~ 6000 stars was selected from an LMC Color 
Magnitude Diagram (CMD) produced from UKST photo- 
graphic plates measured on the APM facility (see Figure]]]). 
These stars were randomly selected with magnitudes and 
colours in the ranges 16 < R < 18 and —1.0 < Bj — R < 2.5, 
in order to have as few biases as possible on the sample se- 
lection. The fields were chosen centered on each previously 
microlensed star. 

The present dataset is the result of 3 observing cam- 
paigns from 1999-2002, using the 2dF (2 Degree Field) in- 
strument at the 3.9m Anglo- Australian Telescope. This in- 
strument is a fibre spectrograph coupled to a robot posi- 
tioner that is capable of observing the spectra of up to 
400 sources simultaneously over a 2-degree diameter circular 
field. Due to adverse weather conditions only 5 field setups 
were obtained, all in UK Schmidt Telescope field 56, for a 
total of 1576 target stars observed. Figure |2| shows the dis- 
tribution of these sources over the face of the LMC. Table 1 
provides a summary of the observing log. The 1200V grating 
was used, giving spectra with dispersion of 1.1 A/pixel over 
the wavelength range 4625A< A < 5765A. For each fibre 
setup, three 1800 sec exposures were combined. 

The spectral images were debiased, flat-fielded and 
wavelength calibrated using the excellent 2dF pipeline soft- 
ware. The spectrum of the night sky was monitored with 




Figure 1. The left-hand panel shows a Hess diagram of the LMC 
field. The right-hand panel shows the portion of CMD where 
our sample stars are selected, from a box of 16 < R < 18 
and —1.0 < -Bj — R < 2.5. The targets are marked accord- 
ing to the radial velocity measured in this survey: stars with 
-100 km s -1 < v h < lOOkms -1 are plotted as filled circles, stars 
with 100 kms~' < vy t < 170 kms" 1 displayed with "star" sym- 
bols, stars with 170kms -1 < < 380 kms - 1 are plotted as 
dots. The two remaining stars with > 400 kms" 1 are marked 
as open circles. 




Figure 2. The location of the target objects in the LMC. The 
stars are marked according to their radial velocities, as in Fig- 
ure 1. 



^ 20 fibres in each setup, and mean sky spectrum was sub- 
tracted from the object spectra. Judging from the residuals 
at the position of the bright [OI]5577A sky line, the sky 
subtraction process is accurate to better than ~ 1%. After 
combining the 3 exposures, the average quality of the spec- 
tra of our survey stars was greater than S/N = 10. The 
wavelength calibration was checked by fitting the peak of 
the [OI]5577A sky line. We estimate that the typical veloc- 
ity error is ~ 15 ± 10 kms" 1 . 

Due to the wide colour range of the stars observed 
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Figure 3. The heliocentric radial velocity distribution of stars in 
the sample. 



Table 1. Summary of Observations 



Date 


Field 


RA 


DEC 


Epoch 


Nov 24 2000 


F056 Conf 05 


05 04 03 


-69 33 18 


2000 


Jan 5 2002 


F056 Conf 01 


05 14 44 


-68 48 01 


2000 


Jan 5 2002 


F056 Conf 02 


05 17 14 


-70 46 58 


2000 


Jan 5 2002 


F056 Conf 04 


05 26 14 


-70 21 14 


2000 


Jan 5 2002 


F056 Conf 21 


05 24 03 


-68 49 12 


2000 



( — 1 < B,t — R < 2.5), the cost of observing a sufficient num- 
ber of radiaf vefocity standard stars of spectraf types cover- 
ing this cofour range woufd have been prohibitive with 2dF. 
Instead we chose to use survey stars as radial velocity stan- 
dards, picking 10 high S/N stars at roughly even intervals 
of Bj — R colour over the above range. Each survey star 
was crosscorrelated against each template star, giving a rel- 
ative velocity measurement, the corresponding un certainty, 
and Tonry- Davis cross-correlation "R" parameter (Tonry & 



Davis 
"R' : 



I97£) 



The template spectrum that gave the highest 
value was considered to provide the best crosscorrela- 
tion match with a particular survey star, and the resulting 
radial velocity difference was used to compute the heliocen- 
tric radial velocity of the survey star. The zero-point velocity 
of one of the (K-type) LMC radial velocity "standards" was 
measured with respect to a genuine radial velocity standard 
of the same spectral type. The radial velocity offsets of the 
remaining 9 LMC radial velocity "standards" were calcu- 
lated by fitting the peak of the numerous LMC population 
(where we assume that there is no systematic offset in ve- 
locity between stars of different spectral types in the LMC). 

To ensure a good quality sample, we applied three pa- 
rameter cuts to the dataset: we require the Tonry-Davis 
cross-correlation parameter to be R > 5; the heliocentric 
radial velocities are constrained to have \vh\ < 1000 km s -1 ; 
and the radial velocity uncertainties (determined from the fit 
of the cross-correlation peak) were set to be Svh < 50 km s~ 1 . 
The resulting dataset, cleaned of low signal to noise spectra 
and spectra that were not well matched by our library of 
10 "standard" stars, comprises 1347 objects. The velocity 
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Figure 4. Heliocentric radial velocity as a function of star colour. 
Upper panel: data. Lower Panel: predictions of the Galaxy model 
with a color conversion Bj — R = 1.25(B — V). Note the outliers 
in the data cannot be explained by the Galaxy model. 



histogram of this sample is displayed in Figure q. A scatter 
diagram is also shown in Figure W for the velocity and the 
color of our sample stars. 



3 THE EXPECTED GALACTIC AND LMC 
POPULATIONS 

The LMC is located at £ = 280°, i.e. roughly toward the di- 
rection of anti-rotation, and is below the disk at b = —33°. 
Since the Sun rotates about the Galactic Cent re 10kms~ 
faster than the LSR (Dehnen & Binney 1998), most local 



stars will be seen towards lower heliocentric radial veloci- 
ties. However, the main effect is that the projection of the 
circular velocity vector of the disk decreases beyond the tan- 
gent point. Thus the majority of Galactic disk stars seen 
along the line of sight towards the LMC will be seen with 
Vh < Okms -1 . In Figures ^(lower-panel) and ^ we display a 
kinematic star-counts model, which contains disk, thick disk 
and spheroid populations. The expected number of sources 
per square degree within our magnitude range is 1012 (disk: 
dashe d), 620 (thi ck disk: dot-dashed) and 370 (spheroid: dot- 
ted) jTbatej |l994| ). 

Note that, of the stars with Vh < 100 kms , most have 
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Figure 5. The expected velocity distribution (upper panel) and 
color distribution (lower panel' ) of Ga l actic stars, according to the 



[bata (1994). The velocity contri- 



Galactic kinematics model of 
butions of the disk, thick disk and stellar spheroid are displayed, 
respectively, by a dashed, a dot-dashed and a dotted line. The 
solid line shows the sum of the three components. 



intermediate colours B — V ~ 1, and are therefore typically 
Solar-luminosity main-seq uence stars with Mr ~ 6 (e.g. 

The distance range correspond- 



Binney & Merrifield 



199S 



ing to our magnitude selection range 16 < R < 18, is then 
~ 1 kpc < d < 2.5 kpc. Given the Galactic latitude of the 
fields, these disk stars are located several scale heights be- 
low the disk, and mostly beyond the tangent point. Nearby, 
intrinsically fainter, and redder disk dwarfs hardly make it 
into the magnitude selection region. 

The simplest explanation of our data is as follows: 

(i) The 40 stars with Vh < 100 km s" 1 are normal disk 
dwarfs, and display an expected radial velocity distribution. 

(ii) The 1291 stars with 170 km s" 1 < v h < 380 km s" 1 
are drawn from a single simple Gaussian distribution, with 
intrinsic dispersion a v = 24kms~ . 

(iii) The 3 stars with 100 kms" 1 < v h < 170 kms" 1 and 
Bj — R < have suspect velocities because their radial ve- 
locities as derived from different (but similar) templates do 
not agree. 

(iv) The 11 stars with 100 kms" 1 < v < 170 kms -1 and 
Bj — R > are spheroid stars. 



(v) The 2 stars with v > 380 kms 1 are halo stars. 

According to the model displayed in Figure ^, the 
number of stars in the velocity range 100 kms" 1 < Vh < 
170 kms" 1 should number only 5% of the number in the ve- 
locity range —100 kms" 1 < Vh < 100 kms -1 . The observed 
number is 35%. Is the excess of ~ 11 stars due to an extra- 
neous unexpected population that could account for the ob- 
served microlensing? This population, however, would only 
represent a small fraction of the mass required to account 
for the observed frequency microlensing. Hence this data ef- 
fectively rules out the existence of an extraneous kinematic 
population that exceeds 1% of the LMC. 



4 COMPARISON WITH PREVIOUS WORK 

The distribution of radial velocities of LMC stars has been 
studied by many authors. The kinematical properties of the 



ing HI (Kim et al. |l998 


e.g.), star clusters (Freeman et al. 


198S ; Schommer et al. 


1992), planetary nebulae (|Mcathcr- 


ingham et alj 1988), HII regions and supergiants (Feitzinger 


et al. 1977). Recently van der Marel et al (2002) have com- 



bined the carbon stars dataset of Kunkel et al. (1997). which 
covers the perip hery of the LMC at a ll position angles, with 

J ^OOJli w hich were selected 



Hardy_ et al 



ED 



573 stars from 

fro m the surveys of Blanco et al. (1980). van der Marel et 



(2002) describe the radial velocity field with a rotation 
of stars in the LMC centered on a C M = 5 h 27.6 m ± 3.9 m 
and 5cm = -69° 52' ± 25', which is about 1.2° ± 0.6° away 
from the gas kinematical center. The kinematics of stars are 
roughly described by a flat rotation of 50 kms -1 outside 4° 
with an inclination i = 34.7° ± 6.2° and a position angle 
= 129.9° ± 6°. The velocity dispersion for the carbon 
stars is 20-22 kms -1 between 1 and 3.5 kpc from the center, 
followed by a decline to 16-17 kms -1 between 3.5 and 7 kpc 
from the center, and a subsequent increase to 21-22 kms" 1 
bet ween 7 and 9 kpc. The velocity dispersion increases with 



age ( Gyuk et al 



2000 



e.g.), with a range from a ~ 6kms _1 
for the youngest populations (e.g., supergiants, HII regions, 
HI gas) to a m 30 kms" 1 for the oldest populations (e.g., 
old long-period variables, old clusters). The carbon stars are 
part of the intermediate- age population which is believed to 
be fairly representative for the bulk of the mass in the LMC. 
Our samp le appear to be cons i stent with the carbon star 
sample of 



van der Marel et al 



(|00§). Figure I is made to 



examine the question whether there is any significant veloc- 
ity variation with spatial position. Given the sparseness of 
our sample, we do not see any strong systematic variation 
across our fields. Our fields are also mostly on the short-axis 
of the LMC disk, hence the signature due to rotation is also 
very weak. This supports our treatment of our sample as a 
whole. 



5 ANALYSIS 

We now investigate the question of whether we could have 
detected a significant population superimposed on the LMC. 

Figure W is a diagram similar to the Hess diagram, but 
shows the Galactic model projected in the radial velocity 
and colour plane. We normalise the Galactic model using 
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Figure 6. A scatter plot of the observed velocity vs. the po- 
sition angle (from the North and centered on the LMC center). 
Large symbols denote velocity outliers. The two dashed sinusoidal 
curves show the projected rotation pattern of the LMC disk (pre- 
dicted from models of van der Marel et al. 2002 at r = 1.5° from 
the LMC center), but offset by 4cr where a ~ 24kms _1 is the 
velocity dispersion of our LMC sample). Note that our sample 
fields do not cover all PA's uniformly, and our fields are mostly 
on the short-axis (PA= 129.9°) of the LMC disk. 



those neighbouring field histograms. Compare this with the 
distribution of possible Galactic stars in our sample e.g. the 
4 blue stars with 50 < v < 150 km/s. It is difficult for a 
standard Galactic model to account for the velocity distri- 
bution we are observing even for the redder objects in the 
< v < 150 km/s range. 

We also test the Galactic model by comparing it with 
a colour histogram from our plate data in that region in the 
16 < R < 18 magnitude range. With plausible assumptions 
about completeness, Figure [j] shows the color histograms 
in a few UKST fields #54, #56 and #58. Except for the 
reference field #56, which is on the LMC, these fields are 
sufficiently away from the LMC to represent the Galactic 
distribution in the general direction of the LMC. While there 
are significant differences between #54 and #58, the overall 
distribution resembles the prediction from our Galaxy model 
(cf. Fig. ^), i.e., most of Galactic stars in the LMC direc- 
tion have a largely bimodal color distribution peaked around 
Bj - R = 1.25(B - V) = 1 and 1.5. Our Galaxy model, al- 
beit simple, should not significantly bias our conclusion. 



6 CONCLUSIONS 

If the microlensing in the LMC is indeed due to self-lensing, 
then kinematic outliers are expected towards the LMC in 
their distributions of radial velocity, proper motion, pro- 
jected distance from the center of the LMC, distance mod- 
ulus and reddening (Zhao 1999a,b). The lack of features in 






Figure 7. From top to bottom color histograms for UKST fields 
#54, #56, and #58. Field #58 is toward the Galactic Plane at 
about (l,b) = (281°, -25°) and #54 at (286°, -42°) with E(B — 
V) = 0.11 and 0.04 respectively. The middle panel is for the 
reference field #56, which is roughly centered on the LMC at 
(281°, -34°) with E(B - V) = 0.23. We used a colour cut which 
removes the stars in the blue and late M in the red. 



our radial velocity sample argues tentatively for a relatively 
smooth and uniform distribution of stars in the LMC with 
very few outliers. Our data does not support the notion that 
there is a large fluffy stellar halo around the LMC (e.g., 
Evans & Kerins 2000). 

On the other hand, it is very unlikely for the hypothet- 
ical (Milky Way) halo white dwarfs to account for the lens- 
mg in the LMC. |Zhac| ( ^002| ) argues that any violent galac- 
tic winds following early epoch of star bursts would signifi- 
cantly weaken the potential wells of galaxies. The Milky Way 
galaxy would have been disintegrated if more than half of its 
dynamical mass were blown off violently in a wind generated 
by the formation of numerous macho progenitors. Since the 
Milky Way galaxy is very tightly bound, and there is no sign 
of any globular cluster escaping from the Milky Way either, 
these should imply an upper limit on the baryons participat- 
ing the early star bursts and baryons locked in stellar rem- 
nants, such as white dwarfs. The white dwarfs should not 
make up more than 1-5% of the total mass of the Galaxy. 
Similar arguments also imply upper limits for the amount of 
neutron stars and stellar black holes, in galaxy halos. This 
dynamical upper limit is exceeded by the amount of halo 
white dwarfs claimed in recent proper motion searches and 
in earlier microlensing observations (Alcock et al. 2000) in 
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the Ga lactic halo, suggesting these interpretations of obser- 



Olsen K.A.G . & Salyk C. 2002, AJ in press (astro- 



vations to be problematic theoretically. 



We propose that the most likely one of existing models 
for microl ensing in the LMC is perhaps the "mutual- lensing" 
models of Zhao & Evans (2001), where the LMC stars are on 
two distinct planes with different inclinations or warped Q so 
that stars on the front plane could lense those in the plane 
~ lkpc behind. One of the planes may be the LMC disk, 
and the other may well be the irregular bar of the LMC, 
which is known to be offset from the LMC disk. Such a un- 
virialised configuration could only be a transient phase of 
the general process of merging, after perhaps recent interac- 
tions with SMC and the Milky Way. In this type of models 
we do not expect the stars in the two planes to share the 
same rotation around the LMC disk, but the differences in 
systematic velocities should be small because they share a 
common potential well. Our data can exclude any system- 
atic offset of > 70kms _1 safely, but not yet exclude the 
possibility of smaller velocity distortion due to different in- 
clinations in the LMC. A larger sample with more accurate 
velocity would be needed to reveal small differences among 
disk or bar stars of the LMC, and set stronger limits on 
star-star lensing in the LMC. 

We thank Brian Boyle for undertaking the service ob- 
serving at the A AT. 
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1 The LMC disk is known to have a poorly denned inclination 
observationally, and is perhaps warped (e.g., Olsen & Salyk 2002) 



